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Abstract — The steady heat conduction through a circular 
cylinder shell is studied theoretically by applying the method 
of variable separation. The cylinder shell is under conditions 
of constant properties, constant inner wall temperature, 
variable outer wall temperature both in length and azimuth. 
Expressions for temperature distribution and heat flow are 
developed. The characteristics is discussed under different 
variations of outer wall temperature. 
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1. Introduction 

Circular tubes have long been widely used in many 
industries of heat and mass transfer because of its high 
structural strength, easy handling and convenient operation. 
Up to now there have been extensive studies on heat transfer 
through circular tubes in serious conditions such as periodic' 11 , 
layered materials 12 ' 31 , convective coupled 141 , or innovative 
means' 5 ' 61 . Circular tubes have two chief objectives in heat 
transfer engineering. One is to serve as heat transfer section 
surfaces to realize heat energy exchange between two fluid 
medium flowing inside and outside the tubes respectively. In 
this case, the main task is to enhance the heat transfer capacity 
of the tube section as much as possible. The other objective of 
tubes is to be apparatus of transporting some wanted fluid 
medium from one place to another. In this condition, it is 
desired to eliminate leakage of the medium out and minimize 
heat transfer between the tube and the surroundings in order to 
maintain the medium properties constant. 

In the process of heat transfer through circular tubes, 
temperature in the tubes varies generally not only along the 
tube length but also along the radius and the azimuth as well. 
So the heat transfer occurring in tubes in almost all cases are 
of problems of unsteady and three dimensional. However, 
many cases could be simplified to be steady and two or even 
one dimensional problems according to concrete situations. 
Indeed, the steady heat transfer models neglecting temperature 
variations in the tube azimuth and length have simple and 
straight forward solutions, and have been written down in 
many heat transfer books' 7 ' 8 ' 91 long before. In many situations 
however, such negligence might causes severe loss. For 
instance, in the heat transfer of fluid flowing perpendicular to 
a circular tube, the temperature of the inner tube wall at each 
cross section might be considered to be invariant in azimuth. 
Yet for the outer tube wall, the temperature might changes 
remarkably not only along the tube length, but also obvious 
along the azimuth as well because of the heat transfer between 
the outer wall and the outside surroundings. In this situation, 
above mentioned models in present text or handbooks could 



result in large errors in predicting the tube wall temperature 
and in designing new innovative compact heat exchangers' 101 . 

In this study, a theoretical investigation is carried out on 
the heat conduction through a circular cylinder shell in steady 
state in which the temperature varies only along the tube 
length on the inner tube wall but changes both along the tube 
azimuth and length on the outer tube wall. The heat 
conduction in the tube length or axial direction resulted from 
the temperature gradient in the tube length is taken to be the 
effect of an inner heat source. By using variable separation, 
expressions for temperature distribution in the tube shell and 
heat conduction rate through the tube are obtained. Effects of 
the temperature variations in azimuth of the outer wall are 
discussed. 

2. Formulation 

2.1 Problem 

Consider a problem of steady heat conduction through a 
circular cylinder shell shown in figure 1 with Cartesian and 
polar coordinate systems as well. The inner and outer radius of 
the tube wall are r\ and r 2 respectively. The tube has a length 
of L and a constant heat conductivity of L The temperature of 
an arbitrary point r,8, z within the shell is denoted as t(r,6 ,z). 
The inner wall temperature keeps constant in azimuth but 
varies along the tube length and is denoted as t x (z). The outer 
wall temperature is symmetry about x coordinate, varying 
along the length and azimuth and is denoted to be t 2 {6 ,z). The 
temperature of the left and right boundary walls are denoted to 
be h(r,8) and t^(r,9) respectively. Under such conditions, we 
treat the heat transfer in tube length to be resulted from an 
inner heat source of q v in the tube wall. Then the governing 
energy equation with its boundary conditions are as follows. 

1 9 * d 2 t d 2 t q v 

(r — ) + + + — = ( 

r dr dr r 2 % 2 dz 2 * 

1) 

r = r x , t=h(z) (2) 
r = r 2 , t=t 2 (0,z) (3) 
z =0, t=t 3 (r,0) (4) 
Z=L, t=t 4 (r,0) (5) 
dt(r,e=<d,z)ld6=dt(r,6=n,z)ld6. (6) 

2.2 Simplification 

2.2.1 Treatment of Temperature 

As mentioned above, the temperature of the inner wall of 
the tube is assumed being constant along azimuth but varies 
along length. Still more, it is further assumed that the axial 
temperature 
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Fig. 1 Circular Cylinder Heat Conduction 



variation of the inner wall represents meanwhile the 
temperature variation of every point in the cross section at the 
same axial location as the inner wall. Under such conditions, 
the temperature at point r, 9, z within the tube shell could be 
written as 

t{r,e,z)=h{z)+T{r,0). (7) 
Substituting Eq.(7) into Eqs.(l),(2) and (3) yields 

1 



d ar 

(r — ) + — 

r dr dr - 2 



d 2 T 



d 2 t 



r 39 dz 

r=n, r=o 



X 







(8) 
(9) 
(10) 

(11) 



r=r 2 , T=t 2 {9,z)-h(z)=m . 
By applying Eq.(7), Eq.(5) - Eq.(4) gives 

U{r,9)-t 3 (rM=tm-tm ■ 
Equation (10) implies that t\(z) and h(9,z) have same change 
in axial direction. Equation (11) indicates that with 
negligence of temperature difference in both radius and 
azimuth at the tube cross section, the temperature change 
between left and right boundary surfaces could be determined 
by the temperature change along the tube length. 
2.2.2 Inner heat source 

The component of heat conduction through the tube length 
is considered to be the result of an inner heat source, by which 
the overall heat generation rate is 



/br(r 2 2 



hm) ^{r 7 ?-ff)(fi(L)-fi(0)) 



(12) 

Where f 3m and f 4m denote the mean temperature of the left and 
right boundary walls respectively. Then the heat generation rate 
of the inner heat source of unit volume of the solid tube shell 
becomes 

(*!<!,) -*!(())) 



(13) 



3. Solution 



3.1Temperature t x 



Recalling the definition of Eq. (7), decompose Eq.(8) 
and let 



d 2 t x 







(14) 



With the conditions z=Q,h(z)=h(0) and z=L,t\(z)=h(L), the 
solution to Eq. (14) could be 



f 1 (z) = f 1 (0)+ 



fiW-fiCO) , 1 q 



z+- 



<U-z z ) 



(15) 



L 2 X 

3.2 Temperature T 

With the method of variable separation, set 

T(r,8)=R(r)0(0). (16) 
Substituting Eq. (16) into Eq. (8) with accounting for Eq. 
(14), one obtains 



3 2 R 



dr 



1 St d 2 
.0+ + -^ R = Ot 



divided by Eq. (16), Eq. (17) is deduced to 



2 a 2 B 

r o R 



r OK 
R dr 



d 2 



(17) 



(18) 



Rdr* R dr <p^ 

Where n is an arbitrary constant number. After some 
rearrangement, two ordinary differential equations are 
deduced to be 

/2, 



2 d^R dR 2d n 
r — — + r n R — 



dr 1 



dr 



d 2 
'dO 2 



(19) 
(20) 



The boundary conditions to Eq. (20) are 

d@ d0 

= 0, = 9 = n, = and 0(9+2k rr) = 0(9) , 

d9 d9 

(21) 

where k is an arbitrary integer number. The third expression 
in Eq. (21) means that what ever the solutions might be, it 
would be the same point and thereafter the temperature 
would be the same. 

Deducing from Eqs. (20) and (21) with n=0, one gets 
00= c o (22) 

where c ois an unknown constant. When n=l,2,3 , we 

have 
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n =c On cos(n0). (23) 
Solving Eq. (19) with «=0 yields 

Rq =c l0 +c u lnr . (24) 

When n= 1 ,2,3 , there gives 

R n =c 2n r n +c 3n r- n . (25) 

From boundary condition Eq. (9), one get R=Q at r=r { . 

Substituting it with some rearrangement, Eqs. (24) and (25) 
become respectively as 

Rq =Ciq— Cyolnr/lnri (26) 



2n—n 



R n = c 2n r " ~ c 2n r l " r 



(27) 



then 



T(r,ff) = Y/n(r,&)= Tfln(r)& n (0) 



n=0 



= Cq(1-Iii r / In r x ) + ^c :l (r n —r^ n r n )cos(n9), (28) 



n=l 



where c =CooCio, c n =Co a C2n- All of the constants c are arbitrary and 
need to be evaluated with the other boundary condition Eq. (10). 
After substitution Eq. (10) into Eq. (28), one arrives at 



T(r 2 ,0) = c (1 - In r 2 / In r x ) + 2c n (r 2 n - rf n r 2 n ) cos(n 9) = f{9). 

n=l 

(29) 

= J_r_ f(0) 

7TJ0 1 



Then 



c 



-lnr 2 l lnr x 



-de 



(30) 



2 f- f{0) 



_ 2 r 



2n — n 

r, - r, r- 



cos(n0)d0 n=l,2,3. 



(31) 

'2 ~ '1 '2 " 

By applying Eqs. (15) and (28) into Eq. (7), the temperature 
of every point within the circular cylinder shell could be 
calculated. 

3.3 Rate of Heat Conduction 

In present problem, the rate of heat conduction through a 
circular cylinder shell consists of the heat conducted between the 
left and right boundary walls and the heat conducted between the 
inner and outer walls. The former component is dependent upon 
the overall slope of the temperature variation in the axial 
direction and could be evaluated by Eq. (12). The latter 
component, i.e. the heat conduction in the radial direction 
through the inner and outer walls is 



Qr 



fit Pit 

-22 — rd0- 
J0 dr 



C 1 dT 
-2/1 — rd0 
JO dr 



■ 7,7th l(cq Anri ) . 



(32) 



4. Discussion 

In practical cases, when the temperature of the inner wall, 
outer wall, left and right walls are given, c and c„ could then be 
calculated with Eqs. (30) and (31). The temperature of every 
point within the circular cylinder tube could be evaluated by 
using Eqs. (13), (15) and (28) into (7). The heat conduction rate 
in the tube length could be calculated by Eq. (12), the heat 
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conduction rate in radius direction could be calculated by 
Eq. (32). When the temperature of the inner and outer walls 
change very slightly along the azimuth so that could be 
ignored, only the change along the radius is apparent, then 
f 2 (#,z) - h(z)=j{0), becomes a constant and is denoted as 
f(9)=At, one gets 

c =- — ; — — — and c n =0 (33) 



\ — lnr 2 /ln r x 
Q r = 2nXL- 



At 



(34) 

ln{r x /r 2 ) 

This result reduces to that in general textbooks about 
one-dimensional heat conduction through circular 
cylinder in radial direction. 

If the inner wall temperature keeps invariant while 
the outer wall temperature changes lineally along 
azimuth and expressed as 

t 2 (0,z)- t l (z)=f(0)=At (l +s9. 
Here e is a constant coefficient. Then it is obtained as 
At 1 ST 

c o - 



\-lnr 2 / lnr x 2\-lnr 2 / In r x 
2s (-!)"-! 



(35) 



™2 n 
7m r 2 



2n — n 
h r 2 



Q r = 27tX L 



At + 1/2S7T 



(36) 



(37) 



ln{r x / r 2 ) 

Thirdly, when the outer wall temperature changes 
non-lineally along azimuth while the inner wall 
temperature still being invariant, take the temperature 
difference between the inner wall and outer wall, for 
example written as 

h(0,z)- h{z)=f(9)=At o +p0 2 . 
Here /? is also a constant coefficient. Then we have 



c o ; 



1 



/3tt 2 



1 - lnr 2 I lnr x 3 1 - lnr 2 I lnr x 



Ap (-1)" 



2 

n r 



Q r = 



Ar +l/3^ 2 



(38) 



(39) 



(40) 



ln(r x lr 2 ) 

At these results, one can concludes that when the outer 
wall temperature varies while that of the inner wall keeps 
constant along the tube azimuth, temperature of every 
point within the solid cylinder will definitely differ from 
that of the radial heat conduction with the inner and outer 
wall temperatures both invariant along the azimuth. As 
for the radial heat flux, there also seems apparently to be 
of significant differences. In order to make this question 
clear, further analysis is needed. 

The question to be considered is that, while the 
inner wall temperature is azimuthal constant, the outer 
wall temperature varies lineally or non-lineally in 
azimuth. No matter how the outer wall temperature 
changes, it has an equivalent mean value along the 
azimuth. That is the outer wall temperature varies along 
the azimuth equivalently above or down an same value. 
For convenience, let the difference between the invariant 
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inner wall temperature and the mean temperature of the outer 
wall be At. 

When the outer wall temperature changes linearly, in 
accordance with above presumption, we have such relation 



\X& Q +e6)&6= [\t&0 . 

JO JO 



And then 

At +l/2s x =At. (41) 
When the outer wall temperature changes quadratically in 
azimuth, we have 



^(At + /?<9 2 )d<9 = ^Atd0 . 



Then it is got 

At a +l/3fi x 2 =At. (42) 
Recalling Eqs. (37) and (40), we see that it is obvious that 
the rate of radial heat conduction remains the same for the outer 
wall temperature changing in the azimuth linearly or 
quadratically as long as their mean temperature being the same. 

From this point, in the conditions of invariant inner wall 
temperature along azimuth, we may deduce that no matter how 
the outer wall temperature varies along the azimuth, so long as 
the mean temperature keeps the same, the radial heat flux 
conducted though a circular tube remains the same. As for the 
temperature within the tube shell, present approach could play 
an appropriate role. cq,c\,c 2 ,...c„ should be calculated by Eqs. 
(30) and (31) with related inner and outer wall temperature 
substituted, the number of constants should be chosen according 
to practical situation. After all, temperature could be evaluated 
by using Eqs. (7), (13),(15) and (28) of this study. 

Under present treatment, heat conduction in tube length and 
tube radius don't interfere each other, but exert additional effects 
in temperature. 

5. Conclusion 

In this study the steady heat conduction in a circular 
cylinder both through radial and axial directions is investigated 
theoretically for the conditions in which the temperature of the 
inner wall keeps constant but that of the outer wall varies along 
the azimuth, apart from this, both the inner and outer wall 
temperature vary equivalently in axial direction. The method of 
variable separation is used in which the heat conduction in axial 
direction is treated as the effect of an inner heat source. Taking 
the inner wall temperature as starting point, the partial energy 
equation is separated into two sets of ordinary equations, and 
solved respectively. The result illustrates that under present 
conditions, the temperature distribution within the circular 
cylinder is determined by the variation of the outer wall 



be calculated by the means provided in present text or 
handbooks with the mean outer wall temperature used in 
the temperature difference. 
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temperature both in azimuthal and axial direction. It is 
indicated that when the outer wall temperature changes along the 
azimuth around a mean value, the radial heat conduction could 
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